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Abstract 

In the attempt to derive the regression theorem from the fluctuation dis- 
sipation theorem several authors claim the violation of the former theorem 
in the quantum case. Here we pose the question: does it exists a quantum 
fluctuation dissipation theorem (QFDT) in its conventional interpretation ? 
It is shown that the relation usually called as the QFDT is the condition 
of detailed macroscopic energetic balance. Following this interpretation the 

existing conflict between the two theorems in the quantum case is removed. 
PACS numbers: 03.65. Ca, 05.30.-d, 05.40.+j 
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I. INTRODUCTION 



Under thermal equilibrium conditions the behavior of fluctuations of macroscopic ob- 
servables of a physical system is governed by relationships which are formulated usually in 
terms of the regression theorem (the so called Onsager hypothesis [1]) and the fluctuation- 
dissipation theorem [2-6] (also known as the Nyquist relation). The former pertains to the 
time domain and states that the relaxation of a correlation of fluctuations is described by 
the same law governing the irreversible processes of the observable quantity itself. The lat- 
ter pertains to the frequency domain and interrelates in some universal way the spectral 
characteristics of fluctuations and hnear response (i.e. dissipation) of an observable of the 
physical system. Often, the fluctuation dissipation theorem is written as: 



where S^xij^) is the spectral density of fluctuations of the observable x, Im{ax{uj)} the 
imaginary part of the generalized susceptibility which is responsible for the relaxation, and 



is the Planck relating factor describing the average energy of the thermal radiation of the 
fleld mode with frequency ui. 

In classical case /i — > 0, the two theorems give the same description of the spectrum of 
thermal fluctuations. By contrast, within the more general quantum case there appears a 
conflict between these two theorems since they predict different behavior of Sxxi}^)- The 
quantum regression theorem (QRT) claims that the spectrum of fluctuations is determined 
only by the permitted transitions between energetic states of the system. In turn, the 
quantum fluctuation dissipation theorem (QFDT) states, that besides the eigenfrequencies 
of the system (as it required by the QRT) in the spectrum of the fluctuations there exist 
additionally the so called Matsubara frequencies, Q„ = i^^^n, where n — ±1, ±2, .... The 
origin of these frequencies is related with the poles of the Planck factor g{uj). 



(1) 



g^{u)^hcoth{^)^2n{n + ]^) 



(2) 
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This conflict is usually interpreted as a violation of the QRT (see, for example, Refs. 
[7-9]). In its most evident form such a violation is demonstrated in Ref. [9], where the 
conclusion statement announced that " there is no quantum regression theorem" . The proof 
of the general character of such a statement is based on the fact that the violation of QRT 
follows from QFDT. However, a proof that QRT is vahd independently of QFDT was given 
by Lax [10] on the basis of the general principles of quantum statistics (see also Refs. [11- 
13]). Since in Refs. [8,9] it is claimed that QFDT contradicts the vahdity of QRT, we argue 
that the origin of such a conflict is related with QFDT and its interpretation (see also Ref. 
[14]). 

The aim of this article is to address this issue by considering the origin of such a conflict 
from a formal mathematical point of view. 



II. FORMAL DERIVATION OF THE QFDT 

In the operator representation, the symmetrized correlation function of the fluctuations 
of an observable x is written as [4-6]: 



CUr) = -Tr{pJx(T)x(0) +x(0)x(t)]} 
while the corresponding linear response function is given by Kubo formula [4-6]: 



(3) 



a.(T) = -e(T)Tr{p,[x(T)x(0) -x(0)x(t)]} 



(4) 



where ps is the density operator which describes some stationary state of the physical system 
under test characterized by Hamiltonian Hgi and ©(r) is the unit step function. 

The linear response described by Eq. (4) implies that the interaction between the system 
and the radiation is determined by the semi-classical Hamiltonian V — —xf{t), where f{t) 
is a classical force. By performing the Fourier transform of Eqs. (3) and (4), in accordance 
with ref. [6] one obtains: 





1 - 


.hl'm{ax{u!)} . 


~ 2 . 



(5) 



where 

/oo 
Tr{psx{T)x{0)}e"^'^ = 27r ^ - uj) (6) 

'°° m,n 

is the spectral density corresponding to the one-directional in time correlation function, 
Tr{psx{T)x{0)} and Jxx{—^) corresponds to Tr{psx{0)x{T)}, pn is the probability to find 
the system in the eigenstate with energy E^, x^n is the matrix representation of the operator 
X and cumn — {Em — En)/h is the frequency associated with the transition between the 
energetic states Em and En- 

Under thermal equilibrium ps = exp{—(5Hs)/Tr{exp{—(3Hs)} with (5 = l/{kT). For the 
derivation of the QFDT an explicit expression for Jxx{'^) is not necessary, it is sufficient the 
existence of the quantum spectral relation [6] : 

Jxx{-^) = e-^'^Jxxiuj) (7) 
By substituting Eq. (7) into Eq. (5) one obtains: 

Jxx{^^) (8) 
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.fiIm{ax{oj)} . 


~ 2 


l-p{uj). 



where p{ijj) — exp{—l3hw). According to Ref. [3-6], the QFDT is then obtained by excluding 
the factor Jxxij^) which is common to both Sxxi}^) and Im{ax{uj)} ■ As a consequence of 
such a derivation, the zeros of Im{ax{uj)} determined by the factor 1 — p{ijj) — (see Eq. 
(8)) become the poles of the Planck factor in Eq. (1), i.e., they originate the Matsubara 
frequencies and, in turn, the QRT-QFDT conflict. Thus, the poles of the Planck factor at the 
Matsubara frequencies can not be considered as independent of the frequency dependence 
of Im{ax{oj)} which, in accordance with Eq. (8), has zeros at the same frequencies: 



Im{ax{uj)] 



= (9) 

U!=Q.n 



Therefore, from a mathematical point of view, in this case the right-hand side of Eq. (1) 
contains an indefinite form of ^-type. 
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Now we pose the following open question: does the constraint given by Eq. (9) merely 
represent the formal requirement necessary for a rigorous derivation of the QFDT, or bring 
a proper physical meaning ? 

An attempt to answer this question is detailed in the following section. 

III. INTERRELATION BETWEEN THE QFDT AND THE PRINCIPLES OF 

ENERGY BALANCE 

As physical model we shall consider a sufficiently large isolated system subdivided into 
two subsystems. The first corresponds to some physical system under test, the second 
represents the surrounding world. In this case the total Hamiltonian can be written as: 
H — Hs + Ht + V, where Hs and Ht are Hamiltonians of subsystems, and V describes the 
interaction between these subsystems. By using the standard procedure [15] to construct the 
master equations for the statistical operators of each subsystem, {i = S,T), and assuming 
that the interaction is weak one obtains the following equation for time variations of the 
average energy < Hi > in the i-th subsystem: 



1 roo ^ ^ ^ 

n Jo 



where < >- TrdpA}, V{-t) = exp[-^iHs + HT)r]Vexp[^iHs + Ht)t]. 

The energy exchange described by Eq. (10) satisfies the conservation law for the total 
energy of an isolated system in the form corresponding to the assumption of a weak interac- 
tion: j-^ < Hs + Ht >— 0. By assuming that the interaction Hamiltonian V — —xf is hnear 
and factorized with respect to the variables of both subsystems the matrix representation 
in the right-hand side of Eq. (10) takes the form: 



d 
dt 



<Hs> 
< Ht> 



S T 



m,n M,N 



mn 
T 



{PmPM - PnPN)^{^mn + ^Mn) (H) 



Prom Eq. (11) one directly obtains the condition of the microscopic detailed energetic 
balance (MiDEB): 
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pipl=pX or 4 = ^ (12) 

Pn Pm 

which must be satisfied with respect to only those energy states of both subsystems that are 
directly involved in the interaction, i.e., when lu^^ — ujjjj^, ^ 0, I/matP 7^ 0. 

Now, let us formulate the conditions of the energy balance at the macroscopic level of 
description. For this sake in Eq. (11) wc replace the term 5(t<;^„j + ^mn) by / 5{uj^^ — 
u!)S{uj]^j^ + u!)du! and rewrite Eq. (11) by using the matrix representation of the asymmetric 
spectral density J{u;) given by Eq. (6): 

^ < >^ < Ht >^ ^ J uj[J^^{-uj)Jff{uj) - J^^{uj)Jff{-uj)\duj (13) 

Prom Eq. (13) we obtain the condition of macroscopic detailed energy balance (MaDEB) 
as: 

Jxx{-uj)Jff{uj) = Jxx{^)Jsf{-^) (14) 

which requires to be fulfilled for any value of the current frequency ui. We notice that the 
condition given by Eq. (14) is not the only form which can be used to express such a detailed 
balance. For example, by using the definitions of Im{ax{oj)} and S^xi'^) given by Eq. (5) 
it is easy to show that Eq. (14) can be rewritten in an equivalent form as: 

Yf{u)Im{ax{u)} = Y,{u)Im{af{u)} (15) 

where for Yi{ijj) one can use any of the three spectral densities, namely: Jn{co), Jn{—Lj) and 
Sii{u!) (i = X, f). The equivalence of Eqs. (14) and (15) from the point of view of the formal 
determination of J(a;), S{u!) and Im{a{u!)} given by Eqs. (5) and (6) does not mean their 
entire macroscopic equivalence. This is related to the fact, that with respect to Eq. (14), 
Eq. (15) involves additionally the notion of the macroscopic function of linear response, 
ai{u!), all properties of which cannot be obtained from Eqs. (5) and (6) only. 

Let us consider the conditions when the MiDEB given by Eq. (12) and the MaDEB given 
by Eqs. (14) and (15) can be considered as equivalent descriptions of the energy balance 
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under steady state. In the following we shall use the formulation of balance conditions 
through some ratio of characteristics of a subsystem [see the right-hand side version of Eq. 
(12)] since such ratios sometimes are universal functions which are independent from the 
internal properties of the interacting subsystems. Let us rewrite Eq. (14) in the form: 



= p{uj) or Jii{-uj) ^p{uj)Jii{uj) (16) 



where p{u!) is common for both subsystems factor which is some single-valued function of 
the current frequency and it satisfies the condition p{—oj) — p~^{u!). It is easy to show that 
the equivalence of the energy balance description given by Eq. (12) and (16) is satisfied if 
p{u!) can be defined from the microscopic level as: 



Pn Pm 



(17) 



Let us rewrite the MaDEB condition given by Eq. (15) in a form analogous to Eq. (16) 

= ^^(a;) or r,(u;) = 5^(u;)/m{a,(u;)} (18) 



Im{ax{uj)} Im{af{Lj)} 

Here the function of current frequency g^{uj) will depend on which spectral density [i.e., 

the symmetric S{u!) or asymmetric J(±a;)] is used to formulate the balance conditions. 

The fulfillment of Eq. (16) allows us to represent the frequency dependence of Su^uj) and 

Im{ai{u!)} of both subsystems in a form entirely analogous to thermal equilibrium [see Eq. 

(8)] with the only difference that now p{u!) is not necessarily given by the thermal value. 

From the above it is easy to see that all the functions g^{(^) in Eq. (18) are determined by 

by the frequency dependence of p{u!) only: 

'1 , /or J{u;) 

g^^iu) = piu) , for J{-u) (19) 

1 - p{uj) 

^+p{u;) , for S{u;) 
When p{u)) is a universal function of frequency [e.g., in thermal equilibrium when p{u!) — 
exp{— f3hu)] there is the possibility to formulate the MaDEB conditions in terms of ex- 
pressions which relate in some universal way the macro-characteristics of only one of the 



interacting subsystem [see the right-hand side expressions in Eqs. (16) and (18)]. However, 
in so doing, it is easy to loose the physical meaning of these expressions. 

In thermal equilibrium, Eq. (18) allows the QFDT to be given an alternative physical 
interpretation with respect to the conventional one. Indeed, by replacing in Eq. (18) Yi{u!) 
with the symmetrical spectral density Sii{u) one obtains the MaDEB condition [right-hand 
side expression in Eq. (18)] in a form which is identical for both the subsystems and which 
coincides with the conventional form of the QFDT given by Eq. (1). This allows us to 
conclude that the QFDT describes the detailed energetic balance between the interacting 
physical systems under thermal equilibrium. As a consequence, the usual interpretation that 
the frequency dependences of the Planck factor g^{u;) [see Eq. (2)] and of the imaginary part 
of the generalized susceptibility Im{ai{uj)} are independent (the source of the QRT-QFDT 
conflict) is in contradiction with both the MiDEB and MaDEB principles. When these 
principles are fulfilled, the poles of g^{u!) and the zeros of Im{ai{u!)} are determined by the 
same factor 1 —p{uj) [see Eqs. (19) and (8), respectively]. The neglect of this property when 
treating fluctuation phenomena corresponds in essence to the violation of the MiDEB and 
MaDEB principles and, hence, it implies the violation of the condition of stationarity at 
least at the Matsubara frequencies a; = Q„ where p{fln) — 1- 

The formulation of the balance conditions represented in terms of a ratio of macro- 
characteristics of one or of both subsystems can serve as a source of incorrect interpretation 
of the MaDEB principle and, as a consequence, can lead to the violation of the energy conser- 
vation law for subsystems with partially overlapping energetic spectra. The typical example 
is the interaction between the harmonic oscillator with eigenfrequency cUg and the thermal 
bath characterized by a quasi-continuous spectrum. By using for the matrix elements of 
oscillator the relations [15]: |a;j„nP = (l/'^s)7s ^(^ + l)<^m,n+i and uj^^ = uJsirn — n), from 
Eq. (6) one obtains asymmetric spectral densities J(±a;) represented as: 



J{u) 



7^(0;) h 



N{uj) + 1 



Sius - u) (20) 



where N{u!s) — S^o^Pn ~ [exp{Phu!s) — 1] ^, 7^(<^s) is the electro-dipole matrix element 



square which can depend upon the oscillator frequency. In Eq.(20), due to the properties 
of the 5-function we have performed the replacement usually used in such situations when 
dependence of 7^(<^s) and N{u!s) on cug is replaced by their dependence on the current 
frequency cu. In this way, from Eq. (20) one obtains: 

J{uj) ' N{u) + l ^ ' 

which defines in explicit form p{uj) in the whole frequency range, while the spectrum of a 
single oscillator is defined at the frequency ui = uis only, i.e. it is described by a 5-function 
which in Eq. (21) disappeared from explicit consideration. 

By considering the thermal bath, also in this case it is used the harmonic oscillator model 
which describes one mode of the radiation field. To obtain the spectral densities J(±a;) for 
the thermal bath it is sufficient to perform in the right-hand side of Eq. (20) a summation 
over all the possible modes, that is over subindex s. In so doing, the dependence of ^'^{oJs) 
and N{ijjs) on ujs but not on uj can remain. By using the standard procedure for systems 
with a quasi-continuous spectrum, i.e., by replacing by / G{yJs)(i^s-i where Giujg) is the 
density of the radiation mode with frequency ujg in the interval duJsi one obtains: 



J(-a;) 
J{u) 



N{u) ■ 
N{uj) + 1. 



(22) 



Prom Eq. (22) one can again obtain Eq. (21), which is now vahd in the whole frequency 
range. When a single oscillator interacts with the thermal bath, the interaction takes place 
not in the whole frequency range, as it would follow from Eqs. (21) or (18), but only at the 
oscillator eigen-frequency, where the energy exchange is only possible. 



IV. CONCLUSIONS AND OPEN QUESTIONS 

Here we have shown that the QFDT expressed in its usual form of Eq. (1) in essence rep- 
resents the macroscopic principle of detailed energy balance between a physical system and 
the surrounding world interacting with it. Such an interpretation, and the related restrictions 
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concerning the frequency dependence of the imaginary part of the generahzed susceptibihty 
of the system [see Eq. (9)], leads to two main consequences. First, the announced conflict be- 
tween QRT and QFDT related with the Matsubara frequencies is reconciled. Second, some 
conventional notions as macroscopic dissipation and relaxation closely related to the QFDT 
should be revisited. Indeed, these notions cannot be treated as internal characteristics of 
one of the two interacting subsystems which compose the total isolated system. 

In the framework of the macroscopic approach where the energy dissipation is described 
by Im{a{uj)}, the energy conservation law under the energy balance implies that the power 
dissipated by one of subsystems from another must be returned back, i.e. it is equal to the 
power dissipated by another subsystem from the first. The form of the MaDEB given by Eq. 
(15) just express this mathematically, by claiming that the spectra of the power dissipated by 
each subsystem from another are equal. Therefore, with respect to one subsystem, one side 
of Eq. (15) can be treated as the dissipated power and another side as the power returned 
back, i.e. as emitted power. This is illustrated by the following example, which considers 
the system interaction with thermal bath. By using Eq. (22), the MaDEB conditions given 
by Eq. (15) take the form: 



2 ?iuj 



7VM + |±2 



Im{ax{uj)} — OJ 



(23) 



The left-hand side of Eq. (23) describes the power dissipated by the system when this power 
is provided by the thermal bath. The right-hand side of Eq. (23) describes the power which 
the system returns back to the thermal bath. If the only photon part huuN^uj), of the full 
energy of the field is involved into the absorption process, then the returning spectrum, 
Jxxi—i^), corresponds to the spectrum of the spontaneous emission processes occurring in 
the system. This is easy shown by using the explicit form of J{—(^) given by Eq. (6). If now 
one reverses time, i.e. replace cu by —cu, the energy of the thermal bath which is involved 
into dissipation process will be hu;[N(u;) + 1]. This means that with respect to one of the 
interacting subsystems the processes of energy exchange in cquihbrium conditions shows 
an asymmetry under time reversal. If the system has absorbed the full energy of the field 
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together with the zero-field value, huj\N{uj) + |] (the case of the QFDT), to fulfill MaDEB 
it would be necessary that the system returns back the same full energy. In this case, 
the emitted spectrum which returns back would be described by the symmetric spectral 
density S^xi^^)- Here the following question arises, which spectrum of fiuctuations and 
correspondingly which correlation function are experimentally measured: the asymmetric or 
symmetric ones ? 

In applications, the QFDT is often used to describe the relaxation phenomena in systems 
interacting with a thermal bath. In this interpretation as one of the formulation of the 
MaDEB conditions there appears a question about the correctness of such a use. Indeed, 
if the relaxation process implies that the system is approaching the stationary state, it 
means that the system is not in a stationary state. It is evident that under nonstationary 
conditions all the relations considered above loose their meaning, and any attempt to use 
them will imply a violation of the energy conservation law. Thus, the QFDT can be apphed 
for the description of only those relaxation phenomena which occur under energy equilibrium 
conditions, as described by the QFDT itself. 

The MaDEB conditions lead to another important property of the relaxation processes: 
the description of the relaxation in the interacting subsystems which compose an isolated 
system cannot be considered in independent way. From the definition of the linear response 
function ai{u;) it follows that the spectral densities of the fiuctuations in each of the sub- 
systems must be interrelated as Sii{u!) — \ai{u!)\'^Sjj{u!) where i ^ j. Together with the 
MaDEB in the form given by Eq. (18) this leads to additional restrictions put on the 
frequency dependence of /m{Q;j(a;)} pertaining to each of the two subsystems: 

Im{a;\u;)} = Im{aj{u;)} (24) 

where Im{a^^{u)} — Im{ai{uj)} /\Im{ai{u)}\'^ is the relaxation law in one of the sub- 
systems, which in accordance with Eq. (24) is determined by the imaginary part of the 
generalized susceptibility of the other subsystem. For example, in the presence of an inter- 
action of some systems with a thermal bath characterized by the radiation mode density 
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G{uj) — iiP"K ^ [15] the relaxation in the system must be described by the law: 

Im{a-\u^)} = ^[J//M - = (25) 

If the electrodipole matrix element square, 7^, is independent of a;, one obtains the usual 
law of viscous friction. 
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